
use "$path\datasets\olsdata", clear


local depvar "ch_total_child_places_pc"



ivreg2 `depvar' share_women_offic_data lag_total_child_places_pc , r cluster(gkz )
est store m1
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)


ivreg2 `depvar' share_women_offic_data lag_total_child_places_pc if in_hand_coll_sample==1, r cluster(gkz )
est store m2
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)



ivreg2 `depvar' share_women_offic_data lag_total_child_places_pc if in_estimation_sample==1, r cluster(gkz )
est store m3
sum `depvar' if e(sample)
estadd scalar  mean_depvar =r(mean)
estadd scalar sd_depvar =r(sd)



esttab  m1 m2 m3 using TableA25/tableA25_PanelB.txt, style(tab) replace order( ) mlabel(,none) ///
cells(b(label(coef.) star fmt(%8.3f) ) se(label((z)) par fmt(%6.3f))) ///
collabels(none) ///
  keep (share_women_offic_data lag_total_child_places_pc  ) ///
  stats( N   N_clust mean_depvar sd_depvar , layout( @  @  `""@ (@)""' )  fmt(  %9.0g %9.0g  %9.2f %9.2f  ) ///
 labels(  "N"  "Municipalities" "Mean (SD)"  )) ///
starlevels(* 0.10 ** 0.05 *** 0.01) ///
  varlabels( share_women_offic_data "Share of women"   lag_total_child_places_pc "Child care spots pc(t-1)" )

	








